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This paper discusses the elastic behavior of a very long crosslinked polyelectrolyte chain (Debye- 
Hiickel chain), which is weakly charged. Therefore the response of the crosslinked chain (network) 
on an external constant force / acting on the ends of the chain is considered. A selfconsistent 
variational computation of an effective field theory is employed. It is shown, that the modulus of 
' the polyelectrolyte network has two parts: the first term represents the usual entropy elasticity of 

connected flexible chains and the second term takes into account the electrostatic interaction of the 
monomers. It is proportional to the squared crosslink density and the Debye - screening parameter. 
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,t-» . I- INTRODUCTION 

o 

Polyelectrolytes are of fundamental importance in a wide range of academic sciences to applications. It ranges from 
life sciences such as biology or biochemistry to industrial and practical applications in daily life products. A typical 

S example for the latter are superabsorber materials. These consist of highly crosslinked polyelectrolyte networks which 
are strongly interacting elastic materials. 
>^ | The theoretical interest in polyelectrolytes reaches back to the early days of polymer science (see e.g [Q). Never- 
£h ' theless they belong to the least understood systems in macromolecular science Q , since there are difficulties to apply 
renormalization group theories and scaling ideas in which long ranged (i.e. Coulomb) forces are present. Only very 
, . recently novel types of field theoretic attempts brought progress jjj . 

In the present paper we aim for a theory of the elasticity of polyelectrolyte networks. This is a non trivial task, 
C\| , since most of the classical and modern theories neglect the effect of interactions on elasticity. In neutral networks the 
interactions are mainly given by excluded volume forces. In the dry network state these can be safely neglected for 
most cases, since in such dense systems like polymer melts excluded volume interactions are screened largely 
However, if solvent is added to the network and the network starts to swell, problems arises. Early theories by Flory 
suggested that the elastic part of the free energy and the solvent part, i.e., a Flory Huggins type term, can be added. 
Later on this concept has been named the Flory - Rehner - Hypothesis JjJ in the context of swelling experiments. The 
[ — ' network state would be determined by the minimum of the total free energy. Indeed such approximations are used in 
ON ■ a wide range of application for rubberlike materials in the swollen state || . The comparison with experiments seem 
to be reasonable, i.e., in all cases the modulus is found to be proportional to the crosslink density, although from a 
theoretical point of view the simple addition of the two parts of the free energy must be wrong Q . It must be wrong, 
Q because a complete formulation of the partition function suggests immediately the appearance of cross terms. In a 
subsequent paper , we will analyse the Flory - Rehner - Hypothesis in the context of thermodynamics and its 
application. There it is shown that corresponds to an approximation which neglect fluctuations completely. In the 
present paper this is not the central point and we restrict ourselves to compute the elastic response of a polyelectrolyte 
network including fluctuations on a variational level. 

In neutral networks, however the approximations seem to be not too bad, because the interactions are relatively short 
. , ranged and weak. Moreover the equilibrium swelling degree is then given by the c* - network. When polyelectrolyte 
networks are considered, we cannot expect that the Flory approximation holds. The interactions are long ranged 
and very strong compared to excluded volume interactions. Here a strong interplay of elastic degrees of freedom and 
interactions must be expected. The reason is very simple: The strong interactions change the physical nature and the 
conformation of a charged polymer chain strongly compared to the neutral chain. The state at rest, i.e., a network 
strand without application of an external force does not contain as many degrees of conformation as the equivalent 
neutral chain. Its conformation ranges, depending on ionization and salt content of the solution, from excluded volume 
behavior to a rod like behavior. Thus a more detailed theory is needed to compute the elastic modulus of charged 



*e-mail: wilder @ mpip-mainz . mpg.de 
t e-mail: vilgis@mpip-mainz.mpg.de 



1 



and highly interacting gels. Nevertheless the Flory approximation has been employed also for strongly interacting 
polyelectrolyte gels in bad solvent to study the phase diagram Jll|] . Moreover also recent suggestions |0| have claimed 
an unchanged modulus for polyelectrolyte networks. We will see later, however, the reasons for these statements. 

Most of the classical network theories rely on " single chain models" . This is to say that the elasticity of the network 
can be roughly computed by studying first the elasticity of a single chain. The elastic properties of the entire network 
is then supposed to be given by the partition function of the single chain raised to the power of the number of 
chains. Such ways of computations hold strongly only for weakly interacting systems. Again in neutral networks 
the interactions, i.e., the excluded volume forces, are weak since they are screened, but in polyelectrolyte networks 
these type of assumptions do not seem to be reasonable. Another drawback of these single chain type theories is 
that they ignore the effect of quenched disorder stemming from the crosslinks completely. Although some progress 
for polyelectrolyte networks has been achieved by single chain theories |l3| , which are based on blob pictur es [R~5| , it 
is necessary to study a network model including the quenched disorder produced by the crosslinks (see e.g. [|14|| for a 
review). It has to be kept in mind that so so far only the limits for large and small screenig parameters are studied. 
Large screenig parameters (in terms of a Debye - Hiickel approximation) correspond to a large salt concentration. 
This regime is relatively unintersting, since it corresponds to the good solvent regime. For small screening the chains 
are significantly stretched, and the chains are in the polyelectrolyte regime ]l5| . For networks, however, these two 
limits allone are not sufficient. In strongly crosslinked networks the chain pieces between two crosslinks can easily be 
of the same order as the Debye - Hiickel screening length. Therefore a more elaborated analysis must be carried out. 

In a previous paper | jl6f we had investigated the coupling between elasticity and interactions in a single chain. This 
model calculation was carried out to show that a strong coupling between elastic (conformational) degrees of freedom 
and electrostatic interactions exists. Of course, unlike as in neutral network systems it is not sufficient to consider 
the elasticity of a single chain and generalize the results to a corresponding network of a large number of such chains. 
In the case of polyelectrolyte networks the chains are strongly interacting with each other. Therefore any assumption 
on weakly interacting chains fails. Nevertheless the computational method used in [|l6| has shown to be useful. The 
results presented there ended up in a two regime blob picture. For small forces the de Gennes |l]| electrostatic blob 
model was recovered. In this regime the chain was relatively easy to deform. For larger forces a change in the elasticity 
was predicted. Then the chain entered in a Pincus regime of the prestretched chain. The results have been confirmed 
by simulation and by a self consistent variational principle. Indeed the method developed in []l6f is appropraite to the 
intermediate regime between strong and weak screening. 

In this paper we extend our considerations concerning the single polyelectrolyte chain Jl^ j to a polyelectrolyte 
network. The simplest version of a polymer network was introduced by Deam and Edwards which consists 
of a very (macroscopically) long crosslinked chain. In our case there is no difference between a network made of 
many chains or one long chain, since we assume to be deep in the solid phase, i.e., well beyond the vulcanization 
threshold. As in the case of a single polyelectrolyte chain we are interested in the force-size relationship of the 
network in a good solvent. In contrast to the classical theories the effects of interactions are now taken into account 
explicitly For simplicity we assume a Debye-Hiickel potential for the electrostatic interaction. It is of the form 
V(r) cx 1/r exp(— r/X). This might not be the best choice to reproduce recent simulation data JH| which had shown 
that the Debye-Hiickel potential is not always a good approximation, but the advantage of the potential is that its 
range is controlled by a single parameter, i.e., the Debye screening length A. Indeed little is known about the interplay 
of elasticity (conformation) and strength of the interactions in the theory of network elasticity. The aim of the present 
paper is to learn something about this interplay. Therefore we apply an external force on the ends of the chain. To 
do so, we employ a variational principle and determine the effective propagator of the chain, which allows statements 
about the interplay of conformation and interactions on the elasticity. In the following we will neglect the effect of 
entanglements completely, since we are mainly interested in the contributions of the interactions to elasticity. The 
results given below will then mainly apply for unentangled networks. To produce simple estimates of entanglement 
contributions tube like models can be employed p3. We expect, however, that the effects of entanglements will not 
be very different from those in neutral networks, as long as their number is given. The main problem would then be 
to compute the mean number of entanglements by the presence of the electrostatic interactions. 

The starting point for the computation of the elastic free energy is the Green function of the crosslinked chain 
without an external force, similar to our calculations of the single polyelectrolyte chain Jll|. The force is treated 
through the analytic continuation of the Fourier transformed Greens function to the complex plane. After having 
introduced a field theory the problem is mapped on a Gaussian field theory with a propagator that formally in the 
Fourier space can be written down exactly by making use of the proper self energy. According to the well known 
Feynman variational inequality the sum of the Gaussian free energy and the mean-value of the interacting potential 
has to be minimized with respect to the proper self energy, which is our variational parameter. This leads to a 
non-linear integral equation for the proper self energy, which can be solved approximately. 

The paper is organized as follows. In the next section we present the underlying model and introduce the crosslinks. 
In section III we formulate a field theory and calculate the variational equation for the proper self energy. In section 
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IV this variational equation is solved approximately. These sections will be written out in more detail, since this is 
- to our knowledge - the first time that a strongly interacting network is investigated by this technique. Thus the 
mathematically interested reader might find the main steps of the computation. The results are presented in section 
V. The paper ends with the discussion of the results. 



II. MODEL AND THEORETICAL BACKGROUND 



The starting point for the field theoretic computation is a network formed out of a macroscopically long chain by 
the instantaneous introduction of a sufficiently large number of crosslinks in the liquid phase (see for example fl7| , |l9|| ) . 
We restrict ourselves to a network, which consists of flexible, weakly charged strands. Consequently the Edwards 
model is an appropriate tool to describe the network ^0|. Therefore let us choose as Hamiltonian for the charged 
chain in aqueous solution: 

o i-Ntat /HrA 2 r Ntot Hr 

pH E [r;f] = - 2 / ds - +p d.sf.- 



2l 2 J Q Jo ds 

bz 2 ,JW rN tot exp{-«|r( S )-r(s')|} 

where r(s) represents the chain conformation in three dimensions as a function of the contour variable s, b = 
e 2 1 ' Ar^tQt-sk-QT is the Bjerrum-length, (3 is (fceT) -1 , where k& is the Boltzmann constant and T denotes the abso- 
lute temperature. I is the Kuhn segment length, z is the monomer charge in units of e, eo is the dielectric constant 
and e r the relative dielectric constant. N to t stands for the bare number of monomers on the chain, f is the external 
force and k^ 1 denotes the Debye-Hiickel screening length A. For the introduction of crosslinks we choose the standard 
way suggested by Deam and Edwards fl7|| . We assume for mathematical convenience four functional crosslinks which 
join two arbitrary segments r(si) and r(sj) along the chain. Of course, the value for the free energy then depends on 
the specific choice of the pairs of monomers, but on macroscopic scale only the statistical average on any crosslink 
configuration is of importance. Nevertheless this requires non - Gibbsian statistical mechanics in the sense that the 
crosslink positions represent quenched degrees of freedom. 

The basic problem for the determination of the free energy of the network is the presence of quenched disorder, 
which is contained in the permanent crosslinks. The formation of a crosslink, i.e., the linkage between two arbitrary 
segments r(sj) and r(sj) represents quenched disorder, since the segments are joined together for all times and for all 
thermodynamic situations. Of course, the experimental relevant free energy F depends on the crosslink configuration 
and crosslink realization S. The actual crosslink configuration S is not known in detail, thus the technical difficulty 
is to average the free energy F(S) over all possible crosslink realizations. To do so, it is generally assumed that the 
corresponding distribution V(S) can be determined. 

Let us perform the outline of the idea in more detail. To calculate the free energy F of the network, we have to take 
the statistical average over all crosslink configurations S. This represents the fact, that F is a self averaging quantity. 

F(N tot ,N c ) = ~k B T J dSP(S)lnZ(S). (2) 

Z(S) is the constrained partition function for a network with the crosslink configuration S, N c is the number of 
crosslinks and V(S) is the crosslink distribution function. Since we assume that the crosslinks are instantaneously 
introduced in the liquid phase, V(S) is yielded by the constrained partition function of the liquid phase, which is 
defined in terms of a path integral as: 

Z®(8) = f Vr(s)exp(-f3H E ) J] S[r( Si ) - r( Sj )], (3) 

where denotes that the i-th and the j-th monomer are close to each other in the liquid phase, which means that 
they can form one of the N c crosslinks. Consequently the crosslink distribution function V(S) is simply: 

V(S) = Z(0){S) (4) 

In the following it appears to be reasonable to assume that the so chosen distribution function does not depend on 
the specific deformation of the network. Note that Z(S) differs generally from Z^°\S). Since we are interested in 
deformations of the network, Z(S) is the partition function of the deformed network. 
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To calculate the free energy F (Eq. (g)) explicitly it is convenient to make use of the so called replica trick JL7) 
This, so far, purely mathematical trick relies on the identity 

111 2 = 



dm 
Define 



m— 



F m (N toU N c ) =-fc B Tln /" dSZ(°)(S)Z m (S), (5) 

where m is the replica index. Eq. (^|) shows the origin of the technical term "replica method". By the use of the 
mathematical trick m copies of the system are produced. The free energy F, which is averaged over the disorder of 



the crosslinks, reads 17 



F(N tot ,N c ) = 9F ^ N * 



dm 



(6) 

m=0 



As in the previous paper, the free energy F is calculated by making use of its relation to the corresponding distribution 
functions and Green functions of the corresponding propagator (see Jl6| for the technical details). 

III. FIELD-THEORETICAL FORMULATION 

In the following we will give an outline of the computation of the network elasticity, i.e., our main aim is to compute 
the low deformation modulus of the polyelectrolyte network. Therefore we start from a concentrated polyelectrolyte 
solution, consisting of one macroscopic chain and the appropriate number of counterions to satisfy the condition of 
electro-neutrality. Then the crosslinks are introduced instantaneously by the process described above. The following 
chapter will be very formal, but we think that it is important to do so, since it turned out that non of the methods 
employed for neutral networks can be used in the present context. The main reason for this is that here we do not 
have the option of formulating the problem in terms of one length scale, i.e., the meshsize, but must take into account 
the range of the interaction. 

For the relevant quantities to compute it is necessary to consider the correlation function G(f , iVtot, f , S) of a 
crosslinked polyelectrolyte chain in replica space, where f = (ro,ri, . . . ,r m ) is the replicated 3(m + l)-dimensional 
end-to-end vector of the chain and S the specific crosslink configuration: 



i-ro(Af t ot)=ro i-ri (AT tot )=n rr m (N to t )=r m 

G(f,7V tot ,7V c ,f,S)= / Vr (s) Vr^s) . . . Vr m {s) 

Jr o (0)=0 Jri(0)=0 J Tm (0)=0 



x exp(-/3F E [f,f]) [] [] - ipfo)] (?) 

p=0 

with p the replica index. The force f is only acting on the ends of the final states (replica index 1 to m) of the chain. 
Thus the replicated Hamiltonian He [r, f ] reads 

bz 2 m A, Ntot exp{ _ K | rp(s) _ rp(s0 | } 



2 ^Jo Jo \r p (s)-r p (s')\ 

The important observation is that the replicated Hamiltonian separates in the different replicas. The coupling 
of the replicas comes into play when the average over the distribution V(S) is performed. If we use the standard 
distribution fl7j| the Green function must be computed upon the effective Hamiltonian 

H = H E -z c l ds I ds'Y[S(r p (s)-r p (s')) (9) 

p=0 



JO 
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The latter equation shows the difficulty of theproblem, i.e., all replicas are coupled. Below we choose a different way 
as suggested by Edwards jl7) and Panyukov jl{|. We must employ field theoretic methods ( as also done in |^9[), but 
the treatment of the field theory is very different, since the symmetry of the problem is not of the same nature as in 
the case of neutral networks. 

It is easy to show that the Greens function in the Fourier transformed replica-space depending on a constant force 
f can be calculated by a zero force Greens function. The force can be reintroduced by the analytic continuation of 
the Fourier-space to the complex plane, which means in detail: 



(10) 



G(k, N tot ,N c , f , S) = / df exp{-i(k - z/?f)r}G(f , N tot ,N c , f = 0, S) 

= G(k(°\ k« - 0f, . . . , k< m > - ij3f, JVtot, N c , f = 0, S) 



Here f is the 3(m + l)-dimcnsional force vector (0, f , . . . , f). This is exactly the same mathematical procedure as we 
had already used in the previous calculation concerning the single polyelectrolyte chain (see [|l6|). As a consequence 
of Eq. (|l(]) in the following we neglect the force term in the Hamiltonian and first calculate a zero force correlation 
function. 

The grand canonical correlation function G(k, /io, z c , f) in replicated Fourier space, where k is the wave vector in 
the 3(m + l)-dimensional Fourier transformed replica space, fiQ is the chemical potential of the monomers and z c is 
the fugacity of the crosslinks, can be calculated by the introduction of de Gennes' zero-component field theory (see 
for example fli| ) 

G(k,/i ,z c ,f) = lim / D^i(k)Vi(-k)exp{-/?if$]}. (11) 



H[^\ in Eq. Cjtt 
space reads 1M 



mm)} = 



is the zero force field theoretical Hamiltonian expressed by the n-component field ip, which in Fourier 

I 2 



|^(q)^(-q) + ^-q 2 ^(q)^'(-q) 



■0(qi)^(q 2 )V'(q3)'0(q4)'5(qi + q2 + q3 + tfc) 



qi,q2,q3,q4 



E 

k=0 



qi,q2 



^( qi )^(q 2 )n^(qi + q 2 ) 



(12) 



q3,q4 



,(*) 



(k) 



(k) 

q 3 



where V^ k \q) is the Fourier transform of the Debye-Hiickel potential in the fc-th replica, if) is a n-component vector 
field and J is an abbreviation for the integral notation J & d q/(2n) d with d the dimension of the vector q. In Fourier 

space G(k, /io, z c ) can be written exactly as: 



G(k,/i ,z c ) 



I 2 - 

-k 2 
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S(k, at) 



(13) 



where S(k, z c ) denotes the proper self energy in replica space. 

Since we do not know the exact proper self energy S(k, z c ), we have to calculate it approximately. Therefore we 
now consider an approximate correlation function £(k, z c ) with an approximate proper self-energy M(k, z c ): 



I 2 - 

G(k, no, z c ) = [ fi + -^-k 2 + M (k, zc) 



To proceed with a variational principle we define the Gaussian Hamiltonian TL by 



ip(-k)g- x (%fx ,z c )if(k) 



(14) 



(15) 



5 



The correlation function C?(k, /j,q,z c ) can be calculated within the zero-component field theory of de Gennes (see for 
example |^9| ): 

G(k,n ,z c ) = lim / P^i(k>i(-k)exp{-/3W[$} (16) 

In this notation the well-known Feynman inequality, which can be taken for the calculation of the approximate proper 
self energy M(k, z c ), is given by: 

F<F+(H-H) n (17) 

where 

/ Vi/j exp{-/3H} 

is the mean-value and T the free energy with respect to H. The right hand side of the inequality (|lj) has to be 
minimized with respect to M . T and (H — TL)u can be written in terms of the correlation function G(k, /iq, z c ): 



71 

^=--V I In 



G(q,V,z c ) (19) 



where V is the volume of the replica space. As can be shown easily the second term of the right hand side of inequality 
©is© 



P<H-n> H =-^V^ M(q, Zc m, Zc) ) + ^(rn + l W V ^ ^ " 

J qi,q2 (q x + q 2 j + k l=1 



z, 



c f„2 



^ 7^ + 2n)7 (J_ £(q,M,OJ (20) 

with Vo the volume of a single replica segment and the volume of the whole replica-space. As we want to determine 
the approximate proper self energy M(k, z c ), this function should be the variational parameter. Consequently the 
general minimization condition reads: 

-{T+(H -H) H ) = (21) 



SM(q, 

where 8/SM(q, z c ) denotes the functional derivative with respect to M(q, z c ). After inserting Eqs. ( |i"8| ) and ( |20| ) into 
Eq. ( pl| ) one obtains 

M(k^ 2 ^ 2 ^+W 



q m + |q 2 + M (q, z c) 
f i i 

4 7r / )7 2 / ^22^ 

i q( o, (q(°) +k(°)) 2 + K 2 M + T(q (0) ' k(1 - m) ) 2 + M (q (0) ' k(1 - m) ) 



+ 4:Trbz 2 m 



1 1 



(i) (q« + k«) 2 + K 2 M + KkW,^ 1 ),^ 2 -" 1 )) 2 +Af(k(°),q( 1 ),k( 2 -'™)) 



+ ^(n + 2) ' ' 



2 iq M+fq 2 + M(q,z c )' 

where the short hand notation \S l --- m ) = (k^, . . . , k' m ^) is introduced. This is a non-linear integral equation for 
M(q, z c ), which in the following has to be solved approximately, since the exact solution is unknown. 
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IV. APPROXIMATE SOLUTION FOR THE PROPER SELF-ENERGY 



In analogy to the calculations on the single chain jl6j we restrict ourselves to small external forces applied on the 
ends of the crosslinked chain. Therefore we make the same ansatz for the proper-self energy as in the case of the single 
chain [ fl6| . The only difference is that the proper self energy in this paper is a function depending on the replica-space 
wave vector q: 



Af(q) = a + aiq 2 + O(q 4 ) 
To start with let M r (q) be M(q) - M(0). Then M r (q) is given by 



(23) 



M r (k) = 47r6z 2 (m + 1) 



CO (q(°) +k(°)) 2 



K 2 fl T 



^( q (0) jk (l...m))2 +Mr ( q (0) jk (l...m) 



q( o> (q(0)) 2 + k 2 Ml . + £(q(0))2 + M r (q(0), 0) 



(24) 



where fi r — jJL + ao. Since we assume a replica symmetric solution for the proper-self energy, we consider the second 
derivative of M r (fc) with respect to at vanishing replica-space wave vectors k, which yields a result for a\. The 
details of calculation are exactly the same as for the single chain [fill: 



2bz 2 



ai 



3Z 2 7T^i r K 



1 + 



kI 



(25) 



This result is valid for f3q/ n < 1 or in terms of the force (3f /k < 1. For details see [nj|. It is important to mention 
that the coefficient a\ is independent of the fugacity of the crosslinks z c . Consequently to study the characteristics 
of the network it is necessary to calculate the constant term of the proper-self energy ao- In the following we neglect 
terms of order n, where n is the number of components of the field tp, since we have to take the limit n — > 0: 



M(k) = 4irbz 2 (m+ 1) 



to, (q(°) +k(°)) 2 + K 2 M +^(q(0),k( 1 -™)) 2 + Af(q(°),k( 1 - m )) 



+ z c 



1 



q M 



I 2 -9 



M(q) 



(26) 



Note, that the right hand side of equation ( pq ) diverges even in the limit m — > 0. Since we started with a discrete 
model, we are not allowed to consider infinitesimal small length scales, which means, that we have to introduce a 
cutoff in the integrals of equation (26). An appropriate cutoff is n (see [fL6|). Conse quen tly it is consistent to substitute 
the ansatz ( p3| ) for M(q) in the integrals of Eq. (p6|), since it is valid for |q| < n fl6| . As we intend to calculate ao, 
all k variables in the integrals of Eq. (|2^) vanish, which means that only absolute values of the q variables occur. 
Therefore we introduce spherical coordinates and neglect terms of order 0(m 2 ) in the second integral after having 
transformed to a dimensionless integration variable. Thus we get the following selfconsistent equation for oq: 



a 



2bz 2 (m + I) 



7T 



dg- 



1 



K 2 ^ 



air 



+ 



2z c (p + a Q )^ n +i 



H/6 + ai 

m+»o q 2 (l + 3mln(<7)) 
dq- 



1 



Q 2 



(27) 



2 3(m+l)^|( m +l) r (|( m + !)) (| + ai ) : 

Note, that the integration-variable q is dimensionless. The integrations in Eq. (^7|) can be performed and the result 



is: 



a 



12bz 2 (m+ 1)^/6(^ + 00) 



7r(6/it + 6ao — 6« 2 ai — n 2 l 2 )\J ' [il 2 + 6/iai + ao^ 2 + 6ao&i 

3bz 2 (m + 1)k 
6fi + 6a — Qn 2 ai — k 2 Z 2 



arctan 



k(/ 2 + 6<zi) 



\/6(n + ao)(^ 2 + 6ai 



(28) 
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2z c {/i + a )^ 



2 3(m+l) 7r i(m+l) r (|( m + 1)) (| + ai ) 



i(m+l) 



/i 2 /6- 



1 — 3m + 3m In 




arctan k 




3mi 
H — dilog 



where dilog(a:) is the dilogarithmic function, defined as 

dilog(x) — 



dt 



Mi) 
i - 1 



(29) 



Starting from expression (28) we neglect terms of order 0(m 2 ). As we consider weakly charged networks we only take 
into account terms up to the order of z 2 . Terms of higher order in the charge z are neglected. We are interested in 
the long ranged limit of the Debye-Hiickel potential. So the next step is to make a series expansion with respect to 
small k including terms of order k 3 . Note, that the sequence of the series expansions with respect to the monomer 
charge z and the inverse Debye-Hiickel screening length k is important. For physical reasons it does not make sense 
to consider first the long ranged limit before calculating the weakly charged case. The result of this analysis is: 



a = 



bz 2 (m + 1)k 
M + a o 

Z c 

ir 2 (n + a ) 



K 2 l 2 



2k 2 1 2 



2 12(^ + a ) 
m7 m \n{lni/ ^Jtt) 
T + 2 



2m 1 
T + 6 



ao. 



(30) 



— + - +C(m 



with 7 0.5772157 Euler's constant. Neglecting terms of order Oq in Eq. ( pl| ) leads to a linear equation in ao which 
can be solved. The result for qq is: 



a 



k ( 2bz 2 (m - 
A 4 V t 
2k z c m 
3tt 2 



1) 



bz 2 {m + \) 2K 2 l 2 bz 2 (m+l) 
2 ^ 9^ 
K 2 z c 7m k 2 z c ln(«://Y / 7r)m k 2 z c 



K 2 l 2 bz 2 (m + 1) 
12/i 



(31) 



4tt 2 



2tt 2 



6tt 2 



Since the expansion coefficients a and a\ for the proper self energy Ajf(k) due to the calculation above are known 
approximately it is possible to write down the Greens function formally: 



i 



/i + a + ^-k 2 + aik 2 



(32) 



V. RESULTS 

From Eq. (|32| ) the grand canonical partition function in replica space under the influence of an external constant 
force / acting on the ends of the chain can be calculated: 

H m ( M , Zc , f ) = 0(k<°> , - z/3f , . . . , k(-) - ipi , n, at) | fc=0 



/i + ao — l -^mf3 2 f 2 — aim/3 2 f 2 

Here we reintroduced the force f according to Eq. ( |Io| ) . The monomer chemical potential and the fugacity of cross- 
links which parameterize the grand canonical partition function in Eq. (|33|), should be expressed in terms of the 
parameters N to t and N c . According to Panyukov and Rabin [jnj the expression F m (N tot , N c ) can be calculated by 
the method of steepest descent in the thermodynamic limit Atot, A c — > oo: 

F m {N tot ,N c )/k B T = - ln3(/x, z c ) - A tot /x + A c lnz c (34) 
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Consequently the fugacity z c of cross-links and the chemical potential fi of monomers can be obtained by minimizing 
the right-hand side of Eq. ( j34[) : 

N tot = - dlnEm =°^ Zc) (35) 

6[l 



and 



N c = 0±B^oM (36) 
a In z c 



in the limit of a vanishing replica index m. It can be shown from Eqs. (pq) and (pQ), that: 



TV = 3z c [i 2 K 3 

8k 3 z 2 W 2 tt + 36k^6z 2 tt + 3k 3 z c /i - 18/z 3 7T 2 - 9k^6z 2 tt 2 - 3k 3 6z 2 Ptt 2 



where TV is defined as the crosslink-density TV c /TV tot of the network. Neglecting the charge per monomer z the 
crosslink density reads: 

3z c flK 3 — 187T 2 /i 3 

Since TV is a positive number z c has to be large enough, namely: 

* > ^ (39) 

At this point we make a series expansion of Eq. ( |38| ) with respect to small /i neglecting terms of order /i 2 which is 
possible, because fj, is connected with the crosslink-density TV, that is assumed to be small: 

TV = /i + C>(^ 2 ) (40) 

To calculate z c we make the following ansatz according to inequality ( |39"| ) with a positive y: 

6ttV 



Zr. 



y (4i) 



Now we substitute this ansatz for z c and the result for the chemical potential /i into Eq. ( |37| ) and calculate y. 
Therefore we again make a series expansion with respect to small monomer charges z neglecting terms of order z 4 . 

- 27r 2 TV 2 (967r 2 TV-366TV7r)z 2 2tt 2 TV 2 / - 3bz 2 l 2 ir 2 - 8bz 2 l 2 n\ , . 

TV = ^ — >— + — 3TV + (42) 

y z Kp yn d \ V J 

36z 2 TVtt 2 - 126z 2 TVtt - 3bz 2 l 2 Tr 2 - 8bz 2 l 2 ir , nr 4 . 

+ a + N + + 0(z 4 ) 

3y 

This equation can be solved with respect to y. Again neglecting terms of order z 4 and only taking into account the 
leading term with respect to a small inverse Debye-Hiickel screening length k the result is: 

y= M4 -: )bz2 + 0(z\ K °) (43) 

K 

Now we can write down the result for z c which reads: 

_ 6tt 2 TV 2 3vr(4 - tt)6z 2 

Moreover it can be shown [|l9| , that the conformational free energy of the network is given by 

_,<91n E m (/x,z, 



(44) 



F(N tati N c ,f) = -k B T- 



dm 



(45) 

m=0 
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If the free energy F is known the force-size relationship is the simply calculated by the derivative of F with respect 
to the external force /: 



(R) = _^WW) (46) 

After having inserted the results for z c and ji in the force size relationship, expanded for small charges z neglecting 
terms of order z 4 and higher, since only weakly charged networks are stable, and considered only terms of leading 
order with respect to small k the force-size relationship reads: 

w = (w^ + w) l,,+0 ^ (47) 

This is again the result for the small deformation regime. The result describes a Hookian law for the force extension 
relation and defines the elastic modulus of the polyelectrolyte network. Note further that this result is valid for small 
forces, i.e. /3//k < 1. Therefore the modulus for the small screening and the low deformation regime of the network 
reads: 



2x -1 



which is the central result of our paper. The modulus depends on the density of the crosslink and on the Debye 
screening parameter. Thus both contributions enter in a significant way. Most striking is that part of the modulus 
stemming from the interactions, which depends on the crosslink density N squared. 



VI. DISCUSSION 



In the previous sections we analyzed the force size relationship of a polyelectrolyte network, which was made of a 
very long crosslinked chain. The method presented, which was developed for a single chain [|16| , is in replica-space 
applicable to polyelectrolyte networks. We considered the network in the long ranged limit of the Debye-Hiickel 
potential. 

An important result is that for small forces and weakly charged networks the response to an external force f is 
proportional to f. The modulus G depends on N 2 , where N is the crosslink-density, which is a surprising result in 
contrast to classical considerations on networks, where the modulus is proportional to N p3|. To discuss the result 
for the modulus in more detail let us discuss it in the form 

G- 1 =p(^L + (49) 

The modulus consists of two terms. The first part, Gn oc ksTN, of the modulus is the term corresponding to 
classical rubber elasticity. It is proportional to the temperature and to the crosslink density. This corresponds to 
the usual entropy elasticity of connected flexible chains. The factor 1/6 appears only from the choice of the special 
representation of the network, i.e., one macroscopic chain and has no specific physical meaning. The second part, i.e., 
G\ oc (k^T N 2 k) / (bz 2 ) stems purely from the interactions. It is not entropy elastic, since the Bjerrum length b and 
the Debye screening parameter depends on temperature. The overall temperature dependence is given by Gi oc T 3 / 2 . 
Since both parts have a really distinguished temperature dependence, they can be separated experimentally in a 
clear way. Moreover the strong difference in the crosslink dependence, Gn oc N, and G\ oc N 2 allows also a clear 
experimental separation. 

It is interesting to note that the two different terms combine in as two springs in series, one entropic one (the rubber 
network) and another energetic one, coming from the interactions. If the strength of the springs is very different, 
naturally in such systems always the weaker dominates the main elasticity. 

It is important to realize that interactions and elasticity interplay in a clear way. Our results yield then the 
conclusion that the Flory assumption, i.e., adding the different parts of elastic and interaction part, is no longer valid 
in these systems. These approximations are perhaps on a level of the random phase approximation, but clearly the 
field theoretical variational technique used in the present paper is beyond perturbative methods used so far in neutral 
P,p7[ and polyelectrolyte jl2j networks. The Flory - Rehner assumption normally uses for the total free energy of the 
network F = F c i ast i C + F mt , i.e., the addition of the elastic and the interaction part of the free energy HJtJ. In this 
hypothesis the bare expression of the elastic free energy is used which is proportional to the crosslink density N. Here 
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we have shown that due to crossterms the elastic modulus becomes renormalized by the interactions. Thus we claim 
that the simple addition theorem is no longer valid. 

The next step will be to apply the results from the present paper to several experimental situations. In a subsequent 
and less detailed paper fiof we will study the effects on swelling and the thermodynamic behavior. In following works 
we will also study the free energy functional at smaller length scales which will provide informations on the scattering 
behavior of the network. This is especially important when deformation processes on different scales are considered. 
Extensions to entangled systems will also become important within this context. And finally, we have to revisit the 
Debye - Hiickel approximation. So far, we had assumed that the counterion are freely distributed, an assumption 
which corresponds to the Debye - Hiickel approximation j24). Of course, correlation effects will change the picture 
and these higher order effects must be the subject of subsequent studies. 
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